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METHODS AND SYSTEMS FOR ESTIMATING FORMATION 
RESISTIVITY THAT ARE LESS SENSITIVE TO SKIN EFFECTS, 
SHOULDER-BED EFFECTS AND FORMATION DIPS 

Background of Invention 

Field of the Invention 

[0001] The invention relates generally to well logging using resistivity logging 
tools. More particularly, the invention relates to methods and systems for reliable 
determination of formation conductivities. 

Background Art 

[0002] Electromagnetic (EM) induction logging instruments are well known in the 
art. These instruments are used to determine the electrical properties 
(conductivity, or its converse, resistivity) of earth formations penetrated by a 
wellbore. Measurements of formation conductivities may be used to estimate the 
fluid contents of the earth formations because hydrocarbon-bearing earth 
formations are associated with lower conductivity (higher resistivity). 

[0003] Physical principles of EM induction logging are described in, H. G. Doll, 
"Introduction to Induction Logging and Application to Logging of Wells Drilled 
with Oil Based Mud" Journal of Petroleum Technology, vol. 1, p. 148, Society of 
Petroleum Engineers, Richardson Tex. (1949). Since then, many improvements 
and modifications to EM induction logging instruments have since been devised. 
For example, U.S. Patent No. 3,510,757 issued to Huston and assigned to the 
assignee of the present invention discloses an EM logging instrument equipped 
with transverse antennas for measuring formation relative dips. U.S. Patent No. 
6,304,086 Bl issued to Minerbo et al. and assigned to the assignee of the present 
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invention discloses EM logging instruments for evaluating formation resistivities 
in high-contrast thin layers or formations with high dip angles. 

[0004] In a conventional EM logging instrument, the transmitter and receiver coils 
(antennas) have their magnetic dipoles substantially aligned with the longitudinal 
axis of the instrument. The antennas on these tools are referred to as longitudinal 
magnetic dipoles (LMD) antennas. With an LMD tool, eddy currents are induced 
in the earth formations to flow in ground loops that are substantially perpendicular 
to the axis of the tool. The eddy currents then induces secondary magnetic fields 
that in turn generate voltages in the receiver antennas. The magnitudes of the 
detected signals relate to the magnitudes of the eddy currents, which in turn relate 
to the formation conductivities. Certain earth formations, however, consist of thin 
layers of electrically conductive materials interleaved with thin layers of 
substantially non-conductive material In these situations, the responses (signals) 
received by a receiver of an LMD induction instrument will be dominated by the 
eddy currents flowing in the conductive layers. In contrast, less-conductive layers 
will have less or little contribution to the overall responses. Consequently, non- 
conductive layers are often missed using the conventional LMD logging tools 
despite the fact that they are hydrocarbon-bearing. 

[0005] To overcome this problem, new EM induction instruments typically include 

transmitter and/or receiver antennas that have their magnetic dipoles substantially 
perpendicular to the axis of the instrument. These antennas are referred to as 
transverse magnetic dipole (TMD) antennas. The TMD antennas induce eddy 
currents to flow in loops parallel to the tool axis. Thus, the eddy currents flow 
through various layers in a vertical well. In this manner, the sedimentation layers 
act like resistors in series in the eddy current loops. Therefore, the signals 
(voltages) received by a TMD tool is more affected by the more resistive layers - 
i.e., the hydrocarbon-bearing layers. 
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[0006] Examples of TMD instruments include tri-axial induction instruments that 
have three transmitter antennas arranged in orthogonal orientations and three 
receiver antennas oriented in the corresponding orthogonal orientations. In a tri- 
axial induction tool energized in three orthogonal directions, individual receiver 
coils aligned in the same three orthogonal directions measure the voltages induced 
by eddy currents flowing in the surrounding formations. The implementation of 
tri-axial antennas, for example, may be found in U.S. Patent Nos. 3,510,757 issued 
to Huston and 5,781,436 issued to Forgang et al. 

[0007] Tri-axial induction logging instruments can provide improved evaluation of 
heterogeneous rock formations. In addition to being able to locate thin 
hydrocarbon-bearing layers, these tools can also provide improved estimation of 
hydrocarbon reserves in anisotropic reservoirs. Examples of anisotropic reservoirs 
(or formations) include highly-laminated formations. These formations can be 
characterized by two electrical parameters: the resistivity parallel to the bedding 
designated as Rh and the resistivity perpendicular to the bedding designated as R v . 
When wells are drilled in reservoirs that may include anisotropic sedimentation 
layers, it is important for the operators to be able to quickly estimate (preferably in 
real time while the data is being acquired) the degree of anisotropy of a particular 
zone in order to make sure that the well is following the planned path or stay 
within the pay zone. 

[0008] Although tri-axial or TMD instruments can provide improved 
measurements for the evaluation of formation resistivity, the raw measurements 
provided by these instruments, like those obtained with LMD tools, are affected by 
skin effects, environmental effects, shoulder-bed effects, and formation relative 
dips. 

[0009] Skin effects are characterized by non-linear responses of the received 
signals in relation to the formation conductivities. Skin effects result primarily 
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from interactions between eddy currents flowing in adjacent loops in the 
formation. The magnitudes of skin effects depend on a complicated function of 
the coil system operating frequency, the effective length of the antenna system, 
and the conductivity value of the adjacent formation, among other things. 

[0010] Shoulder-bed effects result from eddy currents that flow in sedimentation 
layers lying above and/or below the layer being investigated. The shoulder-bed 
effects are particularly problematic if the layer under investigation is less 
conductive than the adjacent beds. In this case, the conductive adjacent beds will 
have significant (or dominant) contribution in the received signals. 

[0011] To some extent, the skin effects and the shoulder-bed effects can be 
mitigated by tool designs and logging parameters. For example, U.S. Patent Nos. 
2,582,314 issued to Doll and 3,067,383 issued to Tanguy disclose induction tools 
having multiple transmitter and receiver coils arranged in specific relationships to 
"focus" the sonde response function by narrowing the width of the main lobe and 
attenuating the side-lobes. In an alternative approach, U.S. Patent No. 2,790,138 
issued to Poupon discloses an induction logging tool having two separate 
induction coil arrangements, which have the same geometrical center so that 
responses from the two coil arrangements may be used to cancel the contributions 
from the side-lobes. 

[0012] In addition to tool design, signal processing methods have been developed 
to improve measurement accuracy by reducing skin effects and shoulder-bed 
effects. Examples of signal processing approaches include phasor processing 
disclosed in U.S. Patent Nos. 4,513,376 issued to Barber and 4,471,436, issued to 
Schaefer et al. In addition, U.S. Patent Nos. 4,818,946 and 4,513,376 issued to 
Barber disclose methods of processing the induction log measurements to reduce 
the unwanted contributions in the log measurements by minimizing the side-lobes 
in the sonde response function used to correlate the voltage measurements with 
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true formation conductivity. The sonde response function is known as a vertical 
sensitivity curve of the induction tool. Furthermore, U.S. Patent No. 6,304,086 Bl 
issued to Minerbo et al. also discloses a tool and a new processing method (the 
Grimaldi processing) that can provide measurements with minimal skin effects 
and shoulder-bed effects. In addition, this tool can output in real time an estimate 
of R v . and Rh in an anisotropic formation. 

[0013] Although prior art tools and methods can produce good resistivity 
estimates, there is still a need for new methods for formation resistivity evaluation 
that are insensitive to insensitive to skin effects, shoulder-bed effects, and 
formation relative dips. 

Summary of Invention 

[0014] One aspect of the invention relates to methods for determining an electrical 
property of a formation using at least two sets of orthogonal resistivity 
measurements. A method for determining an electrical property of a formation in 
accordance with the invention includes acquiring a first resistivity measurement by 
energizing a first transmitter and receiving a first signal in a first receiver, wherein 
the first transmitter and the first receiver are disposed on the logging tool in a first 
orientation substantially parallel to a longitudinal axis of the logging tool; 
acquiring a second resistivity measurement by energizing a second transmitter and 
receiving a second signal in a second receiver, wherein the second transmitter and 
the second receiver are disposed on the logging tool in a second orientation that is 
substantially orthogonal to the first orientation; and deriving the electrical property 
of the formation from a difference measurement that is derived from the first 
resistivity measurement and the second resistivity measurement 

[0015] One aspect of the invention relates to methods for estimating an anisotropic 
resistivity ratio of an anisotropic formation. A method for estimating an 
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anisotropic resistivity ratio of an anisotropic formation in accordance with the 
invention includes acquiring a first resistivity measurement by energizing a first 
transmitter and receiving a first signal in a first receiver, wherein the first 
transmitter and the first receiver are disposed on the logging tool in a first 
orientation substantially parallel to a longitudinal axis of the logging tool; 
acquiring a second resistivity measurement by energizing a second transmitter and 
receiving a second signal in a second receiver, wherein the second transmitter and 
the second receiver are disposed on the logging tool in a second orientation that is 
substantially orthogonal to the first orientation; and deriving the anisotropic 
resistivity ratio from a ratio of the first resistivity measurement and the second 
resistivity measurement. 

[0016] One aspect of the invention relates to methods for estimating an anisotropic 
resistivity ratio of an anisotropic formation. A method for estimating an 
anisotropic resistivity ratio of an anisotropic formation in accordance with the 
invention includes acquiring a first resistivity measurement by energizing a first 
transmitter and receiving a first signal in a first receiver, wherein the first 
transmitter and the first receiver are disposed on the logging tool in a first 
orientation substantially parallel to a longitudinal axis of the logging tool; 
acquiring a second resistivity measurement by energizing a second transmitter and 
receiving a second signal in a second receiver, wherein the second transmitter and 
the second receiver are disposed on the logging tool in a second orientation that is 
substantially orthogonal to the first orientation; and deriving the anisotropic 
resistivity ratio from a ratio of the first resistivity measurement and the second 
resistivity measurement. 

[0017] One aspect of the invention relates to systems for determining an electrical 
property of a formation. A system for determining an electrical property of a 
formation in accordance with the invention includes a computer having a memory 
storing a program having instructions for: acquiring a first resistivity measurement 
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by energizing a first transmitter and receiving a first signal in a first receiver, 
wherein the first transmitter and the first receiver are disposed on the logging tool 
in a first orientation substantially parallel to a longitudinal axis of the logging tool; 
acquiring a second resistivity measurement by energizing a second transmitter and 
receiving a second signal in a second receiver, wherein the second transmitter and 
the second receiver are disposed on the logging tool in a second orientation that is 
substantially orthogonal to the first orientation; and deriving the electrical property 
of the formation from a difference measurement that is derived from the first 
resistivity measurement and the second resistivity measurement. 

[0018] Other aspects and advantages of the invention will be apparent from the 
following description and the appended claims. 



Brief Description of Drawings 



[0019] 



FIG. 1 shows a prior art well logging system. 



[0020] 



FIG. 2 shows a simple two tri-axial antenna tool that can be used with 



embodiments of the invention. 



[0021] 



FIG. 3 shows the skin effects on various measurements derived from a tri- 



axial measurements according to one embodiment of the invention. 



[0022] 



FIG. 4 shows vertical geometrical factors of various measurements 
according to one embodiment of the invention. 



[0023] 



FIG. 5 shows the radial geometrical factors of various measurements 
according to one embodiment of the invention. 



[0024] 



FIG. 6 shows integrated radial geometrical factors of various measurements 
according to one embodiment of the invention. 
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[0025] FIG. 7 shows formation conductivities derived from various measurements 
in a stair-case formation model having no relative dip according to one 
embodiment of the invention. 

[0026] FIG. 8 shows formation conductivities derived from various measurements 
in a stair-case formation model having a 60-degree relative dip according to one 
embodiment of the invention. 

[0027] FIG. 9 shows formation conductivities derived from various measurements 
in a stair-case formation model having a 90-degree relative dip according to one 
embodiment of the invention. 

[0028] FIG. 10 shows formation conductivities derived from various 
measurements in the Oklahoma 2 formation having no relative dip according to 
one embodiment of the invention. 

[0029] FIG. 11 illustrates better bed boundary definition using a difference 
measurement according to one embodiment of the invention. 

[0030] FIG. 12 shows vertical geometrical factors of various measurements 
according to one embodiment of the invention. 

[0031] FIG. 13 shows formation conductivities derived from various 
measurements in a stair-case formation model having no relative dip according to 
one embodiment of the invention. 

[0032] FIG. 14 shows a cross plot for determining anisotropic ratio of a formation 
according to one embodiment of the invention. 

[0033] FIG. 15 shows a flowchart of a method in accordance with one embodiment 
of the invention. 

[0034] FIG. 16 shows a prior art computer that can be used with embodiments of 
the invention. 
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Detailed Description 

[0035] Embodiments of the present invention relate to methods and systems for 
estimating conductivity of rock formations. Specifically, the present invention is 
directed to methods and systems that can provide consistent and accurate 
conductivity measurements that are less sensitive to skin and shoulder-bed effects, 
regardless of the relative dips of the formation bedding. The conductivities 
obtained using embodiments of the present invention may include an apparent 
conductivity, and a conductivity parallel or perpendicular to the formation 
bedding. In addition, the present invention provides a method for estimating the 
anisotropic resistivity ratio in anisotropic formations. 

[0036] In the following description, the term "conductivity" is used 
interchangeably with "resistivity" because one is the reciprocal of the other and 
either can be used to characterize the electrical properties of a material. In 
addition, the description uses induction logging as an example. However, one of 
ordinary skill in the art would appreciate that embodiments of the invention may 
also be applied to orthogonal sets of measurements obtained from propagation 
logging. Therefore, both induction and propagation measurements are expressly 
within the scope of the invention. Furthermore, methods of the invention may also 
be used to process previously acquired measurements that include two sets of 
orthogonal measurements, i.e., reprocessing of existing data. Thus, any reference 
to logging or data-acquisition steps in this description is only for illustration and is 
not intended to limit every embodiment of the invention. 

[0037] Embodiments of the present invention relate to methods for probing a rock 
formation using at least two orthogonal sets of induction measurements. The at 
least two orthogonal sets of induction measurements are obtained by energizing 
two transmitters having their magnetic dipoles in orthogonal directions (e.g., one 
parallel to and the other perpendicular to the tool axis) and receiving the induced 
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signals (voltages) in two receivers having their magnetic dipoles oriented in the 
same directions as the transmitter magnetic dipoles. For example, the first 
transmitter and receiver pair have their magnetic dipoles aligned with the tool axis 
(i.e., longitudinal magnetic dipole (LMD) antennas), and the second transmitter 
and receiver pair have their magnetic dipoles perpendicular to the tool axis (i.e., 
transverse magnetic dipole (TMD) antennas). One of ordinary skill in the art 
would appreciate that the orthogonal sets of induction measurements may be 
acquired using any induction logging tool equipped with at least one transmitter 
and at lest one receiver, each including an LMD and a TMD antennas. Examples 
of the antenna arrangements may include a tri-axial array, which includes a tri- 
axial transmitter and a tri-axial receiver. In this description, a transmitter or a 
receiver may include a single coil (antenna) or a group of coils arranged in a set. 
For example, a tri-axial transmitter or a tri-axial receiver includes three coils 
arranged in orthogonal directions. While preferred embodiments of the invention 
use orthogonal antennas that are either parallel or perpendicular to the tool axis, 
orthogonal antennas that deviate from these orientations (i.e., tilted antennas) may 
also be used. 

[0038] FIG. 1 shows a schematic of a typical logging system. Certain 
conventional details are omitted in FIG. 1 for clarity of illustration. The logging 
system 200 includes a logging tool 205 adapted to be moveable through a 
borehole. The logging tool 205 is connected to a surface equipment 210 via a 
wireline 215 (or drill string). Although a wireline tool is shown, those skilled in 
the art would appreciate that embodiments of the invention may be implemented 
in wireline or while-drilling (LWD or MWD) operations. The surface equipment 
210 may include a computer. 

[0039] The logging tool 205 may be any conventional induction tool capable of 
providing two orthogonal measurements, e.g., a tri-axial array logging tool that 
includes a tri-axial transmitter and a tri-axial receiver. FIG. 2 shows an exemplary 
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induction tool 220 having a tri-axial array. As shown, the induction tool 220 
includes a tri-axial transmitter 221 and a tri-axial receiver 222. One of ordinary 
skill in the art would appreciate that an antenna may be used as a transmitter or 
receiver. Therefore, specific reference to transmitter or receiver antennas in this 
description is only for clarity of illustration. Furthermore, while a co-located tri- 
axial array, in which the centers of the three antenna coils are co-located, is 
shown, one of ordinary skill in the art would appreciate that other antenna 
configurations may also be used. 

[0040] For clarity, the following description assumes that the tool has a single 
transmitter and a single receiver, as shown in FIG. 2. However, a typical 
resistivity tool may have more than one transmitter and/or more than one receiver. 
A transmitter or a receiver as used herein may include one or more coils arranged 
in a group, such as a tri-axial transmitter or a tri-axial receiver. Furthermore, a set 
of "bucking" coils may be included between each pair of the transmitter and the 
receiver to reduce the mutual couplings between them. The bucking coils 
typically include the same number of coils as that of the receiver coils, and the 
bucking coils are wound in opposite direction to those of the corresponding 
receiver coils.. One of ordinary skill in the art would appreciate that embodiments 
of the invention are not limited by any specific configuration of the resistivity 
logging tool as long as the tool is capable of providing two orthogonal 
measurements of resistivity. 

[0041] In accordance with embodiments of the invention, the first set of 
measurements may be made using the transmitter and receiver antennas (coils) 
having their magnetic dipoles aligned with the tool axis (z-axis). These antennas, 
which are traditionally referred to as longitudinal magnetic dipole (LMD) 
antennas, are referred to as "T Z2 " for the transmitter and "R zz " for the receiver in 
FIG. 2. When T zz is energized, it induces eddy currents to flow in the formations 
in the xy planes, i.e., perpendicular to the tool axis (z-axis). The eddy currents 
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flowing in the xy planes would induce a voltage V zz in the receiver R zz , whose 
magnetic dipole is also aligned with the z-axis. 

[0042] A second set of measurements may be made using transmitter and receiver 
antennas having their magnetic dipoles perpendicular to the tool axis - i.e., along 
the x-axis or y-axis. These transmitter and receiver antennas are traditionally 
referred to as transverse magnetic dipole (TMD) antennas. They are referred to as 
"T xx " or "Tyy" for the transmitter and "R xx " or "Ryy" for the receiver in FIG. 2. 
When transmitter T xx is energized, it induces eddy currents to flow in planes 
perpendicular to the x-axis, e.g., in the yz planes. These eddy currents then induce 
a voltage V xx in the receiver R^, the magnetic dipole of which is aligned with the 
same x-axis direction. One of ordinary skill in the art would appreciate that the 
description for V xx applies equally well to V yy , which is acquired with a 
transmitter and a receiver both having their magnetic dipoles in the y-axis 
direction. Furthermore, any description about or Vyy applies equally well to an 
average (or a weighted average) of and Vyy. 

[0043] The induced voltages (signals) V zz and V^ are proportional to quantities 
shown in equation (1) and (2). 

V zz = K(l - ikL)(l + ikL - -jk 2 ll - ^k 3 L 3 ) 

or V Z2 =K(l + ±k 2 L 2 -U>L 3 +...) (1) 

=iKe ikL (-\ + ikL + k 2 L 2 ) (2 ) 
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where K = —^-r , L is the transmitter-receiver spacing, and k is the wave number. 

In the low frequency limit with negligible dielectric effect, the wave number k 
may be defined as: 

k = JudJio = (3) 

where 8 is the skin effect, which is a function of the operation frequency (a> = 2nJ), 
magnetic permeability of the medium (n), and the conductivity of the medium (o): 



8 = J — — (3a) 
y co fid 

[0044] As can be seen from equation (3a), the skin effect (S) is a function of the 
tool operation frequency (a). If the operating frequency (co) is not too high, then 
the skin effect may not be significant. If the skin effect (S) is not significant 
relative to the transmitter-receiver spacing (L) (i.e., L/S <1), then the skin effect is 
only a perturbation to the measurements. Under these conditions, taking the 
difference between the induced voltages (equation (1)) and V** (equation (2)), 
yields equation (4) or (5): 

V„-V„=KQi + U 2 L 2 +^-k 5 L 5 ) (4) 

or, 

y~ - v *x = 3 - — (a/UTp L 5 + i\o)McrL 2 - — (atjioji L 5 J (5) 
[0045] The real part of equation (4) is proportional to the formation conductivity, 



-Vj=cwGL 2 {\-Uco»oilA (6) 
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[0046] Equation (6) shows that the real part of (V zz -V xx ) is proportional to a with 

an error term proportional to — . Therefore, the real part of (V^ - Vxx) shown 

\8) 

in equation (6) provides a measurement of formation conductivity that is much 
less affected by the skin effect (S). In comparison, V zz , which is conventionally 

used to derive the classical o zl , has an error term proportional to and, 

therefore, is more affected by the skin effect (S). Thus, the real part of (V^ - V^) 
should be able to provide more reliable determinations of the formation 
conductivities. This is illustrated in FIG. 3. 

[0047] FIG. 3 shows correlations between true conductivity of a homogeneous 
medium and the derived conductivities based on various measurements obtained 
using a typical resistivity logging tool equipped with at least one LMD antenna 
and at least one TMD antenna, for example a tri-axial resistivity tool. Curve 31 
shows conductivities derived from measurements; curve 32 shows 

conductivities derived from measurements; curve 33 shows conductivities 
derived from Re(V zz ) - 7w(V xx ); and curve 34 shows conductivities derived from 
Vzz-Vxx. It is known that a transverse measurement (e.g., V^) is more sensitive 
than a longitudinal measurements (e.g., V zz ) to skin effects. This is seen by 
greater deviation from the expected values in curve 31 (V^) than in curve 32 
(V 22 ). 

[0048] Curve 33, which is based on ReiW^yim^V^X has less skin effects than 
either curve 31 (V xx ) or curve 32 (V zz ). This is because some skin effects are 
cancelled out in the difference measurements. It is known that the imaginary part 
of a signal correlates with the skin effect and that V^ suffers more from skin 
effects. Therefore, the ^(V^) term provides a good skin effect correction and, 
thus, the difference measurement 7?^(V zz )-/m(V xx ) is expected to have less skin 
effects. What is surprising is that the difference measurement Vzz-V xx is the least 
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affected by the skin effects (see curve 34), even less than the Reiy^-IrnQJ^) 
measurement (curve 33). 

[0049] In addition to being less affected by skin effects, the difference 
measurement (V zz -V xx ) is also less affected by the shoulder-bed effects. FIG. 4 
shows the vertical geometric factors as a function of z/L (where z is the vertical 
distance from the center of the antennas and L is the transmitter-receiver spacing) 
for the V zz , Vxx and V zz -V xx measurements. Curve 41 is the vertical geometrical 
factor for the V xx measurement. Because measures signals induced by eddy 
currents flowing in yz planes, this measurement is expected to be influenced by 
sedimentation layers lying above and below the vertical location of the antennas, 
hence more shoulder-bed effects. Curve 42 shows the vertical geometrical factor 
for the measurement. Because the V zz measurements detect eddy currents 
flowing in the xy planes, they are expected to have less shoulder-bed effects as 
compared with the V xx measurement. Curve 43 shows the vertical geometrical 
factor for the Vzz-Vxx measurement. Clearly, the V Z z-Vxx measurement (curve 43) 
has less shoulder-bed effects than either the V zz (curve 42) or the V^ (curve 41) 
measurement. The improved geometrical factor is another advantage of the 
difference measurement, V zz -Vxx, as compared with the V^ or the V^ 
measurement. The results shown in FIG. 4 is also consistent with resistivity 
profiles derived from these measurements, shown in FIG. 7-9 to be described later. 

[0050] FIG. 5 shows the radial geometrical factors for the V xx , V zz , and Vzz-Vxx 
measurements as a function of r/a (r is the radial distance into the formation and a 
is the radius of the tool). It is clear that the radial response profile of the difference 
measurement, V^-Vxx (curve 53), is shallower than either the V^ (curve 52) or the 
V xx (curve 51) measurement. A shallow radial response profile is associated with 
a better vertical resolution. Therefore, the difference measurement V zz -Vxx is 
expected to have a better vertical resolution than either the Vzz measurement or the 
V xx measurement, consistent with the results shown in FIG. 4. 
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[0051] FIG. 6 shows the integrated radial geometrical factors for the V^, V zz , and 
V zz -V xx measurements as a function of r/a. The integrated radial geometrical 
factors show that the difference measurement V zz -V xx (curve 63) is shallower than 
the V zz measurement (curve 62). However, the difference measurement Vzz-Vxx 
(curve 63) behaves better (no negative values) than the V xx measurement (curve 
61) in the near wellbore region. 

[0052] The above description shows that the difference measurement V zz -Vxx has 
less skin effects, less shoulder-bed effects, and a better vertical resolution, when 
compared with either the V xx measurement or the V zz measurement. These 
properties of the V zz -V xx measurement predict that the V^-V^ difference 
measurement should be able to provide better estimates of true formation 
resistivities and more precise bed boundaries. 

[0053] FIGs. 7-9 show tool responses based on various measurements to a 
staircase resistivity profile in a formation model. The tool has a 27 in. transmitter- 
receiver spacing. FIG. 7 shows the formation conductivities as derived from 
various measurements in a formation with no relative dip. As shown, the 
difference measurements V zz -V xx (curve 72) produce conductivity estimates that 
closely match the true values (curve 71), i.e., with minimum skin effects and 
shoulder-bed effects. In comparison, the traditional LMD measurements, 
(curve 73), produce estimates that are substantially lower than the true 
conductivities due to the skin effects. In addition, the V^ measurement (curve 73) 
also show more shoulder-bed effects. The skin effects and shoulder-bed effects 
are even more severe in the TMD measurements, V^ or Vyy (curve 74). FIG. 7 
also shows the cross-coupling measurements (curve 75) and V zx (curve 76), 
which do not have measurable signals in this vertical well in a homogeneous 
formation. It is clear from FIG. 7 that the difference measurements Vzz-V xx (curve 
72) produce significantly better estimates than do either the V^ (curve 73) or the 
V xx (curve 74) measurements. 
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[0054] FIG. 8 shows the formation conductivities as derived from various 
measurements in formations with a 60-degree relative dip. As shown, the 
difference measurements V zz -V xx (curve 82) produce conductivity estimates that 
closely match the true values (curve 81), i.e., with minimum skin effects and 
shoulder-bed effects. In addition, the estimated values from V^-V^ 
measurements seem insensitive to relative dips; this is apparent from a comparison 
between curve 72 in FIG. 7 and curve 82 in FIG. 8. Again, the V zz measurements 
(curve 83) are more affected by the skin effects and the shoulder-bed effects. The 
skin effects and shoulder-bed effects are even more severe in the Vxx (curve 84) or 
the Vyy (curve 84') measurements. FIG. 8 shows that the cross-coupling 
measurements (curve 85) and V zx (curve 86) are now measurable at bed 
boundaries. Again, FIG. 8 shows that the difference measurements Vzz-Vxx (curve 
82) produce significantly better estimates (i.e., less skin effects and less shoulder- 
bed effects) than either of the V zz (curve 83), V xx (curve 84), or Vyy (curve 84') 
measurements. 

[0055] FIG. 9 shows the formation conductivities as derived from various 
measurements in formations with a 90-degree relative dip, e.g., a horizontal well. 
The difference measurements V^-V^ (curve 92) produce conductivity estimates 
that closely match the true values (curve 91) even in this highly deviated well. 
This confirms that the V zz -V xx measurements are insensitive to relative dips. In 
contrast, the V zz measurements (curve 93) are influenced by the skin effects and 
the shoulder-bed effects. The skin effects and shoulder-bed effects are even more 
severe in the V^ (curve 94) or V yy (curve 94') measurements. The cross coupling 
measurements V^ (curve 95) and V zx (curve 96) do not produce reliable estimates 
of the formation conductivities. 

[0056] FIGs. 7-9 clearly show that the difference measurements V zz -Vxx produce 
significantly better conductivity/resistivity estimates than do either the V zz or V^ 
measurements. Furthermore, the difference measurements Vzz-Vxx also provide 
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better definition of bed boundaries (or boundary contrasts) because these 
measurements provide sharper resistivity/conductivity changes at bed boundaries. 
The better definition of bed boundaries and better estimates of the resistivity in 
each layer makes it possible to provide better weight for each detected boundary, 
i.e., more accurate weight that is proportional to the contrast between the adjacent 
beds. The better resistivity estimates and boundary definition obtained from V z2 - 
Vxx may serve as inputs for the inversion of more complex data sets, such as those 
obtained with tri-axial arrays. More importantly, the difference measurements 
V zz -V xx can produce reliable estimates regardless of the formation dips. These 
advantages of the Vzz-Vxx measurements have also been observed in real 
formations, e.g., the Oklahoma 2 test formation (FIG. 10). 

[0057] FIG. 10 shows the conductivities derived from various measurements in the 
Oklahoma 2 test formation. As shown the V 2z -Vxx measurement (curve 102) has 
significantly less skin effects and shoulder-bed effects than the V zz measurement 
(curve 103) and the V xx measurement (curve 104). 

[0058] As noted above, in addition to providing more accurate estimates of the 
conductivity, the V zz -V xx measurement also provides better defined bed 
boundaries. FIG. 11 shows that the bed boundaries can be accurate identified 
from the inflection points of the V zz -Vxx curves - e.g., by taking the derivatives of 
the curve. Because the resistivity/conductivity derived from V^-V^ 
measurements in each bed is more accurate and because the V zz -V xx measurements 
provide better definition of bed boundaries, the V^- V xx measurements can provide 
more reliable inputs (initial estimates) for the inversion of full tri-axial induction 
measurements, hence speeding up the inversion process. 

[0059] The above description shows that a simple difference measurement, Vzz- 
Vxx, can provide an improved estimate of formation resistivity with less skin 
effects and shoulder-bed effects. However, FIG. 3 also shows that the vertical 
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geometrical factor of the V zz -Vxx measurements is not perfect. Therefore, further 
improvement should be possible. The geometric factor approach has been used to 
further improve the resistivity estimates according to a general formation: a(pV 2Z 
-Vxx), where a and P are constants. With this approach, it was found that the 
difference function l4(3/2Vzz- V^) provides an optimal vertical geometrical factor. 

[0060] FIG. 12 shows the geometrical factor of the V^^V^-V^) measurements 
(curve 123) as compared with the V zz (curve 121) and V xx (curve 122) 
measurements. It is clear that the shoulder-bed effects essentially disappear in the 
difference measurements Y2(3/2V ZZ -V XK ). 

[0061] FIG. 13 shows the conductivities derived from various measurements in a 
stair-case resistivity model formation. A comparison between the V zz -Vxx 
measurements (curve 131) and the 72(3/2 V^-V^) measurements (curve 132) 
reveals that the V^^V^-V^) measurements have less should-bed effects, but 
more skin effects. Both the V^-V^ measurements (curve 131) and the !/2(3/2V zz - 
V xx ) measurements (curve 132) provide better estimates than the V zz 
measurements (curve 133) or the V xx measurements (curve 134). The same 
phenomena are also observed in a vertical well in the Okalahoma 2 formations 
(data not shown). Thus, either the V^-V^ or the ^^V^-V^) measurement, or 
their variants, can provide better formation resistivity estimates than the 
conventional V zz measurements. The choice between the V zz -V xx measurement 
and the 1 /2(3/2V zz -V xx ) measurement will depend on whether one is more 
concerned with the skin effects or the shoulder-bed effects. 

[0062] Some embodiments of the invention provide more reliable estimates of the 
horizontal and vertical resistivities in anisotropic formations. In a formation with 
a low relative dip, if the logging operation is performed with low frequencies (i.e., 
negligible skin effects), the real part of the V zz measurements is proportional to the 
formation conductivity in the horizontal planes, an, while the real part of the 
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measurements is proportional to the formation conductivity in the vertical 
direction, <r v . Various methods have been reported for deriving the ah and a v 

values from the V zz and V xx measurements. However, the Gh and a v values thus 
derived are not always accurate because the V zz and V xx measurements are 
sensitive to the skin effects and the shoulder-bed effects. In addition, the and 
V xx measurements are sensitive to the relative dips. 

[0063] It is found that the ratio of the real parts of these two measurements, 
V zz /2 V xx , is approximately proportional to the square of the anisotropy coefficient 

(^ 2 ), i.e., 

Therefore, if the horizontal conductivity (a h ) can be estimated, then the vertical 
conductivity (a v ) can be derived from the V ZZ /2V XX ratio and the estimated a h . 

[0064] According to embodiments of the present invention, the V^-V^ 
measurements can be used to derive accurate estimates of a h . The a h thus 
obtained can then be used to derive the a v from the estimated a h and the VJIV^ 
ratio. As noted above, the resistivities derived from the Vzz-Vxx measurements are 
less sensitive to skin effects and shoulder-bed effects. Therefore, the estimates of 
a h and a v thus derived are more reliable. 

[0065] FIG. 14 shows a cross plot of -Re(V xx )/2 J Re(V zz ) versus i?e(V zz ) that can be 
used to estimate the R v /Rh ratio for low relative dip formations. Each curve in the 
plot corresponds to a different R v /Rh ratio. These curves may be obtained from 
simulation. To use this plot to obtain the anisotropy ratio and the Rh, a point 
corresponding to the i*e(V xx )/2i?e(V zz ) and Re(V zz ) values is located on the chart 
(e.g., point A shown in FIG. 14). A vertical line is then drawn from point A to 
estimate the R h value (hence, ah). A horizontal line is drawn from point A to 
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estimate the corresponding R v /R h ratio. Once the R v /Rh ratio and the R h value are 
available, R v can then be determined. 

[0066] Thus, the new measurements Vzz-Vxx and V zz /2Vxx according to 

embodiments of the invention can provide extremely useful information in the 
inversion of full tri-axial induction measurements. These measurements provide 
three important parameters: the precise bed boundary location, an approximate 
initial guess to the inversion algorithm and an initial guess for the anisotropy ratio 
R v /Rh. The more accurate estimates can speed up the inversion process. 

[0067] The above approach to estimating R v /Rh from the and V** 

measurements is feasible only when there is no or minimal dip. Some 
embodiments of the invention relate to method for providing reliable resistivity 
estimates in anisotropic formations even when the formations have significant 
relative dips. 

[0068] Effects of formation anisotropy on resistivity measurements in a 

homogeneous anisotropic formation were first described by Moran and Gianzero. 
See Moran and Ginzero, "Effect of Formation Aniostropy on Resistivity Anisotropy 
measurements" Geophysics, Vol. 44, pp. 1266-1286, (1979). The equations for 
anisotropic formation resistivity calculations derived by Moran and Ginzero are 
referenced to a coordinate system tied to the formation layers. Because resistivity 
measurements may be acquired with tools not perpendicular to the formation 
layers (e.g., a formation with relative dips), these equations are often difficult to 
apply. In U.S. Patent No. 6,584,408 issued to Omeragic ("the Omeragic patent") 
and assigned to the assignee of the present invention, these equations were 
simplified to a reference frame relative to the logging tool. This patent is 
incorporated by reference in its entirely. The Omeragic patent discloses a 
procedure for determining anisotropic formation parameters from tri-axial 
measurements in formations with relative dips. According to one procedure, the 
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horizontal conductivity (a h ) of the formation is first determined from the tri-axial 
measurements by using cross coupling terms. Then, the dip angle (a) is derived 
from the measurements and the estimated cr h . Finally, the vertical conductivity 
(a v ) of the formation is derived from the horizontal resistivity (ah) and the dip (a). 

[0069] The method disclosed in the Omeragic patent first derives the horizontal 
conductivity (a h ) from coupling terms in tri-axial measurements by solving the 
following equation: (t zz -L h )(T^ - T h )=(T a f. As noted above, the horizontal 

conductivity (a h ) of a formation may be more reliably obtained from the V^-V^ 
measurements regardless of the relative dips. Therefore, the horizontal 
conductivity (ah) may be more conveniently obtained from the difference 
measurement (V zz -Vxx). 

[0070] Once the horizontal conductivity (ah) is available, the dip angle (a) can be 
determined according to the following equation that is disclosed in the Omeragic 
patent: 

a = tan- (8) 

L h 2 xx 

where and T x2 are the strike-rotated xx (V xx ) and xz couplings (V^ ), and T h is 

the xx coupling (V xx ) in an isotropic formation having a conductivity q&. The term 
"strike-rotated" means that the reference coordinate system has been rotated to 
remove the strike (azimuthal angle) of the dipping plane. The process of removing 
strike from the measurements is disclosed in the Omeragic patent. An alternative 
to obtaining the relative angle (a) is to use: 

a = tan" 1 Lh ~ T " (9) 
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where T zz is the strike-rotated zz coupling (V zz ), and L h is the zz coupling (V 2Z ) in 
an isotropic formation having a conductivity 

[0071] Equations (8) and (9) may be combined using the formula for the sum of 

tan 1 to give: 



a = 0.5tan — — ^ zjJ ^- 



T h -L h +T m -T m 



(10) 



or, 

2T 2 

a = 0.5 tan" 1 ^ ^ r (11) 

[0072] Either of these two equations may be used to calculate the formation dip 
angles (a). However, equation (11) is more robust in extreme cases (0 or 90°), 
and it is undetermined (0/0) only if there is no anisotropy. 

[0073] As noted above, the horizontal conductivity (ah) may be derived from the 
xz coupling (Vxz) or the V zz -V xx measurement. Note that the horizontal 
conductivity (ah) derived from the V zz -Vxx measurement corresponds to the no 
anisotropy situation (i.e., an isotropic formation where L h = T h ). Therefore, 
equation (11) cannot be used if o h is derived from the V zz -Vxx measurement. In 
this case, the alternative approach is to use following expression for the first guess 
of the relative dip, 



a = tan" 1 A L * (12) 



T -T 



h 



[0074] Once the horizontal conductivity (a h ) and the relative (a) are available, 
these parameters can be used to obtain an estimate of the vertical conductivity a v 
according to the following equation: 



23 



PATENT APPLICATION 
ATTORNEY DOCKET NO. 20.2892; 09428.125001 



4xs 



2 ((2r ft+ zJ-fc + 



f +T'f 

* yy zz // 



(13) 



[0075] 



Equation (13) is simpler than the equation for deriving the a v disclosed in 



the Omeragic patent. However, it should be noted that a is a function of cr v . 
Therefore, the above expression should be run recursively. Furthermore, (2T h +L h ) 
is proportional to crin an isotropic formation. Therefore, equation (13) in essence 
derives anisotropy from the difference between the actual tool reading and what 
the tool would read if there is no anisotropy. The vertical conductivity (a v ) 
estimate obtained from equation (13) may then be used in an iterative solver to 
determine refined horizontal conductivity (<r h ) and vertical conductivity (a v ) from 
a full set of tri-axial resistivity measurements. 

[0076] FIG. 15 shows a flowchart illustrating a method 150 in accordance with one 
embodiment of the invention. Initially, two orthogonal voltage measurements 
(<?.£•> and Vxx) are obtained (Step 151). In a logging operation, these 
measurements will be performed at a series of depths. These measurements are 
typically obtained in the form of a voltage log. One of ordinary skill in the art 
would appreciate that the real (in-phase) parts and the imaginary (quadrature or 
out-of-phase) parts of the signals may be logged separately. 

[0077] The two orthogonal measurements at each depth are then manipulated to 
obtain the desired difference measurement (V^-V^ or !4(3/2V zz -Vxx)) and/or the 
ratio of the measurements (V^V^). (Step 152). The resulting difference 
measurement (V^-V^ or Vi^llV^W^) may be used to estimate the formation 
conductivity (a) and to define the bed boundaries (Step 153). The estimated 
formation conductivity (a) and the bed boundaries may be used as initial inputs 
for the inversion of full tri-axial measurements (step 154). As noted above, 
because the V^-V^ or x M^I2V Z2r W y ^ measurements provide better estimates of 
the formation conductivity and bed boundaries, the inversion of the full tri-axial 
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measurements may be performed more efficiently. Particularly, the estimates 
derived from the difference measurements are not sensitive to formation dips. 
Therefore, embodiments of the invention can provide reliable formation resistivity 
parameters regardless of the formation relative dips. 

[0078] If the formation is anisotropic, the derived formation conductivity (a) 
corresponds to an initial estimate of a formation conductivity parallel to the 
bedding plane (a h ). If the formation dips are not significant, the ratio (yjlW^ 
may be used to provide R v /Rh or the anisotropy coefficient (X) (step 155). The 
anisotropy coefficient (X) together with the estimated ah may then be used to 
derived formation relative dips (8) and/or the formation conductivity 
perpendicular to the bedding planes (a v ) (step 156). Embodiments of the 
invention can provide better estimates of the horizontal conductivity (ah) and 
anisotropy coefficient (X). Therefore, the derived parameters for the anisotropic 
formation are more reliable. 

[0079] If the formation dips are not negligible, an alternative approach to deriving 
formation resistivity parameters is to estimate the horizontal conductivity (a h ) 
using the V^-V^ measurements. Once the horizontal conductivity (a h ) is known, 
the dip angle (a) can then be determined from a full set of tri-axial measurements 
according to methods described above, i.e., equations (8) - (12) (step 155). Once 
the horizontal conductivity (a h ) and the dip angle (a) are available, the vertical 
conductivity (a v ) can then be determined (step 156). 

[0080] Note that the above steps for deriving horizontal conductivity (a h ) may 
further include using the initial estimate derived from the difference measurement 
in an iterative process to solve for more accurate horizontal conductivity (ah) 
using the full set of resistivity measurements (e.g., tri-axial measurements). This 
iterative process, if included, may be performed before the horizontal conductivity 



25 



PATENT APPLICATION 
ATTORNEY DOCKET NO. 20.2892; 09428.125001 

(<r h ) is used together with the R v /Rh ratio or the dip angle (a) to solve for the 
vertical conductivity (a v ). 

[0081] Some embodiments of the invention relate to systems for performing the 
methods described above. A system in accordance with embodiments of the 
invention may be a stand-alone unit for performing methods of the invention or 
may be incorporated into a drilling tool. A system in accordance with the 
invention typically includes a processor and a memory. In some embodiments, a 
system may be implemented on a general-purpose computer having a processor, a 
memory, and may optionally include other hardware. For example, as shown in 
FIG. 16, a typical computer (160) includes a processor (163), a random access 
memory (164), and a storage device (e.g., permanent memory or hard disk) (166). 
The computer (160) may also include input means, such as a keyboard (168) and a 
mouse (161), and output means, such as a monitor (162). Note that the general 
purpose computer is only for illustration and embodiments of the invention may 
take other forms (e.g., integrated in a logging tool). 

[0082] In a system in accordance with the invention, the memory stores a program 
readable by the processor. The program, for example, may include instructions for 
performing the above described methods: obtaining resistivity measurements that 
include at least two orthogonal measurements (for example using a tri-axial tool), 
deriving difference measurements and/or a ratio of the two orthogonal 
measurements, estimating the formation conductivity and bed boundaries, 
estimating the anisotropy coefficient, deriving the horizontal and vertical 
conductivity of an anisotropic formation, and deriving dip angles in formations 
with dipping planes. 

[0083] A system in accordance with the present invention provides new and 
improved techniques to evaluate formation electrical properties, e.g., resistivity (or 
conductivity), bed boundaries, anisotropy coefficient, and relative dips. The 



26 



PATENT APPLICATION 
ATTORNEY DOCKET NO. 20.2892; 09428.125001 

programming may be accomplished through the use of one or more program 
storage devices readable by the computer processor and encoding one or more 
programs of instructions executable by the computer for performing the operations 
described herein. The program storage device may take the form of, for example, 
one or more floppy disks, a CD-ROM or other optical disk, a magnetic tape, a 
read-only memory chip (ROM), and other forms of the kind well known in the art. 
The program of instructions may be in "object code," i.e., in binary form that is 
executable directly by the computer, in "source code" that requires compilation or 
interpretation before execution, or in some intermediate form such as partially 
compiled code. The precise forms of the program storage device and of the 
encoding of instructions are immaterial here. 

[0084] Advantages of the invention may include one or more of the following. 
The methods can provide more accurate estimates of formation conductivities, 
which are less affected by skin effects and shoulder-bed effects. In addition, the 
estimates are not affected by the formation dips. Therefore, reliable results may 
be obtained regardless of the formation dips. Furthermore, the methods of the 
invention also can provide more precise definition of bed boundaries. Thus, 
embodiments of the invention can provide more accurate initial estimates of the 
conductivity and bed boundaries for the inversion of full tri-axial measurements. 

[0085] Some embodiments of the invention provide convenient methods for 
calculating the formation anisotropy coefficient. This coupled with the more 
reliable estimate of the horizontal conductivity of the formation makes it possible 
to derive more accurate vertical conductivity of the formation. In addition, some 
embodiments of the invention provide convenient ways to calculate formation 
relative dips. Thus, even in a formation with relative dips reliable estimates of 
anisotropic formation resistivity parameters may be derived. 
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[0086] While the invention has been described with respect to a limited number of 
embodiments, those skilled in the art, having benefit of this disclosure, will 
appreciate that other embodiments can be devised which do not depart from the 
scope of the invention as disclosed herein. Accordingly, the scope of the 
invention should be limited only by the attached claims. 
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